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Abstract. We consider the massive Klein-Gordon field on the half line with 
and without a Robin boundary potential. The field is coupled at the bound- 
ary to a harmonic oscillator. We solve the system classically and observe the 
existence of classical boundary bound states in some regions of the parame- 
ter space. The system is then quantized, the quantum reflection matrix and 
^ I reflection cross section are calculated. Resonances and Ramsauer-Townsend 

^ 1 effects are observed in the cross section. The pole structure of the reflection 

' matrix is discussed. 

\^ , The study of quantum fields with boundaries has been the subject of much work 

' in recent years. There have been two main methods used to study such systems; 

integrable boundary field theory (see, for example, 0121) and boundary pertur- 
bation theory (see ^Ij and references there in). Both approaches have focused on 
situations where the boundary is a non-dynamic object that does not contain its 
' own internal degrees of freedom. Recently some studies have been conducted on the 

, coupling of both classical and quantum fields to boundaries, or impurities, contain- 

ed ' ing additional degrees of freedom, for instance, dynamic boundary conditions have 

been studied for the integrable sine-Gordon 0]-[7], supersymmetric sine-Gordon [5], 
free fermion fields and several authors have studied the coupling of the massless 
Klein-Gordon field to oscillators, e.g. JHI-CSI- These models can be used in the 
study of such physical systems as excited atoms in cavities and quantum wires con- 
taining impurities. Boundary degrees of freedom are also of interest in the study of 
brane-bulk interactions in braneworld universes, see for example, 14 . 
, In this paper we consider another system consisting of a field with a dynamic 

boundary. This 'toy model' consists of a massive Klein-Gordon field in 1+1 dimen- 
sions restricted to the left half line by a boundary. The field is linearly coupled to 
' a harmonic oscillator at the boundary thus introducing the additional degrees of 

freedom. Although this model is entirely linear, and may therefore appear trivial 
on first inspection, we will see it possesses several interesting, non-obvious features. 
The full analysis of these features, which we are able to undertake due to its lin- 
ear nature, will help in the interpretation of similar effects in other non-integrable 
models and in systems where perturbative methods cannot be applied. 

In section 1 we review previous results for the massive Klein-Gordon field with 
the Robin boundary condition. In section 2 we consider the massive Klein-Gordon 
field on the halfline coupled to an harmonic oscillator at the boundary with a Robin 
boundary potential. The system is solved in the classical regime and boundary 
bound states are observed in some regions of the parameter space. We show that 
the classical solutions of the system can be decomposed into independent modes 
of oscillation. We observe that the existence of boundary bound states combined 
with the requirement that the Hamiltonian be bounded below restricts the possible 
values of the parameters of the theory. We then quantize the field directly from the 
classical solutions expressed as a superposition of the independent modes and find 
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the quantum reflection matrices of the systems from the two point function of the 
field. From the reflection matrix we calculate the reflection cross section, resonances 
and Ramsauer-Townsend effects are observed for some ranges of the parameter 
spaces. We discuss the pole structure of the quantum reflection matrices. At the 
end of section 2 we consider the special case when the Robin boundary potential 
is absent. Section 3 contains a discussion on the main results of this paper, and 
possible directions for future work in the area. 

1. The Massive Klein-Gordon Field with a Robin Boundary 

In this section we consider an example of an exactly solvable boundary field 
theory, the massive Klein-Gordon field on the left half line with a Robin boundary 
potential. This system has been studied previously, see, for example, [T^IH]. 

1.1. Classical System. The Hamiltonian for this system is given by 

J° tf + i {d,<j>{x, tjf + \m^c^{x, tf^ dx + iA0(O, tf , (1) 

where (j){x,t) and T:{x,t) are the field and its conjugate momentum which have 
the Poisson bracket relations t), 7r(y, i)} — 5{x — y) The boundary coupling 

parameter A is assumed to be real. 

From we find Hamilton's equations for the system^, 

dt(P{x,t) = {cj){x,t),H}=7r{x,t) , (2) 

dtTT{,x,t) = {■n(x,t),H} ^dlcj){x,t) -■n?cj){x,t) 

- (9,0(0, t)-fA0(O,t)) . (3) 

From (|21) and Q we find the Robin boundary condition by requiring that 7r(a;,t) 
be continuous at the origin, 

9,</)(0,t) = -A0(O,t) . (4) 

From Hamilton's equations and Q and the boundary condition we recover 
the usual equation of motion for the massive Klein-Gordon field, 

d'^(i)(x, t) = dl(t){x, t) - ■m'^(t>{x, t) . (5) 

The equation of motion lO and boundary condition Q are satisfied by the 'bulk' 
solutions, 

0(x,t) = / I cos(pa;) sin(pa;) 

Jo V P 

a{p) cos{LUpt) + sin(a;pt) j dp , (6) 



where a{p) and b{p) are real functions of p and ujp — \/ m?- + p^ . 

As well as the 'bulk' solutions there can exist square integrable boundary bound 
state solutions. Such solutions can found by allowing the momentum, p, of the bulk 
solutions, given in this case by to take imaginary values, p — ig. Because the 
bulk solutions satisfy the boundary condition (@J and the equation of motion of the 
bulk field (O the imaginary momentum solutions will also satisfy these equations. 
However, such solutions will, generally, cause the field to diverge as x — oo 
and so must be discarded'^. However, for special values of g the divergence of 
the field arising from terms proportional to cosh(£ix) and that arising from terms 



^Here and throughout this paper the notation dx4>{0,t) should be read as the derivative with 
respect to x of tpi^yt) evaluated at x = 0, i.e. a shortened notation for dx4'{^}'t)\x—o- 

■^The field is required to be bounded over its domain, thus solutions which diverge on the left 
half line must be discarded. 
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proportional to sinh(pa;) cancel out on the left half line. At such values of g there 
exists a valid boundary bound state. Such a valid boundary bound state exists for 
the Robin boundary &i g = ~\ provided A is negative^. The boundary bound state 
solution is then 

t) = 6^"= (^g coslujgt) + ^ sin(cjgi) j , (7) 

where ag and bg are real amplitudes and LOg = m? — g^ . Note that there is only 
a boundary bound state for negative values of A, i.e., when the contribution of the 
boundary term to the Hamiltonian Q is negative. Such boundaries are called at- 
tractive, those for which the contribution of the boundary term to the Hamiltonian 
is positive are called repulsive boundaries. Repulsive boundaries cannot support 
boundary bound states of this form. Also note that when g^ > rr? the value of cjg 
becomes imaginary and the boundary bound state solution diverges as i ^ icx). 
The most general solution for the Robin boundary is formed by adding the bulk 
solutions to the boundary bound state solution, if one exists. 

We can now try substituting these solutions into the Hamiltonian (Q. Let us 
first assume that A is positive, i.e. the boundary is repulsive and there exists no 
boundary bound state. Taking each term of the Hamiltonian in turn at i = 0, 
substituting the expressions for the bulk solutions © we find that 

ii=\ -h{pf + -^la{pf dp , (8) 



where 



A2 



: = ^-a(p) + , (9) 

1 12 

Hp) ■■ = ^Kp) (i + ^) • (10) 

When A is negative we must also include the boundary bound state solution 
when we are substituting into the Hamiltonian (Q. This results in an additional 
term appearing in ^ corresponding to the energy contribution of the additional 
mode. The Hamiltonian for the attractive Robin boundary is 



H = I hipf + Icoj-aipr dp + hi + hldl , (11) 



where a{p) and b{p) are given by JHl and (|l()(l and 



b 



" ■ (2|A|)i/2 ' 
b„ 



' ■ (2|A|)i/2 • 

Note that if A^ > ni'^ then is negative. This would mean that the Hamilton- 
ian Ullfl would not be bounded below and thus be unphysical. It is therefore a 
requirement on the parameters of the theory that A^ < for the attractive Robin 
boundary, as the repulsive boundary has no such bound state there does not exist 
an equivalent condition, the Hamiltonian is always bounded below. 



^If A is positive there is no solution which decays on the left half line instead there is a solution 
which exponentially grows as x — > — oo which must be rejected. 
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1.2. Quantum System. As both the Hamiltonian for the repulsive ((SJ and the 
attractive Robin boundaries Hll|l are now written as an infinite sum of harmonic 
oscillators we can quantize the field in the canonical fashion. For each mode we 
have a creation and annihilation operator, 



The 'bulk' creation operators, a^(/9), act on the vacuum to produce bulk particles 
of the field. These particles are located far to the left of the boundary a,t t — —oo 
with momentum p = p they propagate to the boundary from which they reflect. 
At t = oo these particles are once again located far to the left of the boundary 
and have momentum p — —p. The creation operator associated with the boundary 
bound state, a],, acts on the vacuum to create a particle of the field with energy 
lower than that of a stationary bulk particle, this particle is trapped close to the 
boundary at all times. 

The one-particle quantum reflection matrix, or i?-matrix, for the massive Klein- 
Gordon field with a linear boundary can be found directly from the two point 
function of the quantized field. 



The first integral on the right hand side of \12\ is the standard two point function 
of the Klein-Gordon field. This corresponds to a particle propagating between the 
space time points (xi^ti) and (12, ^2)- The second integral corresponds to a particle 
propagating between the same two points via the boundary where it is refiected and 
picks up the momentum dependent amplitude R{p), which is the equivalent to the 
i?- matrix found using any other commonly used definitions, see |16j . The third 
term on the right of (|12|l arises from any boundary bound states of the system, 
such contributions will fall off to zero as xi — > —00 and 2:2 — + — 00 due to the 
localization of the boundary bound state particles close to the boundary. 

Using these creation and annihilation operators, the definition of the quantum 
reflection matrix H12(l and the expressions for the general bulk solutions of the field 
it is simple to calculate the i?-matrix for the Robin boundary to be 



Let us consider the analytic structure of H13(l . There is always a single pole located 
at p = — zA. When A is negative this pole appears in the upper complex half plane 
as shown in figure lA where the black circle indicates the pole and the white circle 
indicates a node. The pole occurs at the same value of imaginary momentum for 
which we have observed a classical boundary bound state. The interpretation of 
this pole is clear, it corresponds to the existence of a quantum boundary bound 
state which is equivalent to the classical boundary bound state. However, when A 
is positive, and the pole appears in the lower complex half plane as shown in figure 
IB, there is no classical boundary state to which the pole can correspond and the 
pole does not correspond to a quantum state of the system either. 





-|- boundary bound state contributions. 



(12) 
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Figure 1 . The analytic structure of the quantum reflection matrix 
for the Robin boundary for (A) A < and (B) A > 0. 



2. Coupling the Massive Klein-Gordon Field to a Boundary 
Oscillator in the Presence of a Robin Potential 

In this section we consider the massive Klein- Gordon field in H-1 dimensions 
restricted to the left half line with a Robin boundary potential. In addition the 
field is linearly coupled to a harmonic oscillator at the boundary thus introducing an 
additional degrees of freedom. This preserves the linear, and thus exactly solvable, 
nature of the theory but introduces new and interesting features into the reflection 
cross section. 

2.1. Classical System. The Hamiltonian for this system is 



H = 



+ ty ] dx 



-A</>(0,t)2 



(14) 



where 7r(a;,i) is the conjugate momentum to the field 4>{x,t), and p{t) is the con- 
jugate momentum to the oscillator q{t). These have the usual equal time Poisson 
bracket relations, {^(x, i), 7r(t/, t)} = 5{x — y), {q{t),p{t)} = 1, all other brackets 
are zero. The parameters of the theory are the mass of the field, m, the Robin 
boundary coupling parameter. A, the boundary oscillator coupling parameter, (3, 
and the natural frequency of the boundary oscillator, /i. In natural units^ these pa- 
rameters have dimensions [m ] = [fi] = [A] = [M], [P] = [M3/2] and for the oscillator 
[q{t)] = [M~-^/2] and [p{t)] — [M^/^]. The parameters m, ^, A and /3 are assumed 
to be real and m and /i are taken to be positive. 

Hamilton's equations for the field and oscillator are 



5t0(a;, t) = t), H} = 7r(x, t) , 

dtiT{x, t) = {7r(a;, t),H} — dl4>{x, t) — m^4>{x, t) 
-5{x) (9,0(0, i) + A(/>(0,t)+/3g(i)) 
dtq{t) = {q{t),H]^p{t) , 

dtp{t) = {p{t).il} = -/30(O,t) - il^q{t) . 



(15) 

(16) 
(17) 
(18) 



From (|16|l we find the boundary condition necessary for 7r(x, t) to be continuous to 
be 



5,0(0, t) = ~A0(O,t)-/3g(O 



(19) 



^Setting c = h = 1, gives relations between the dimensions, [AI] = [L ^] = [T 
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Combining the boundary condition (|19l) and Hamilton's equations (|15() and (|16|l we 
recover the equation of motion for the Klein-Gordon field 

dt(l){x, t) = dl(j){x, t) - m^(j){x, t) . (20) 

Similarly we find the equation of motion for the boundary oscillator from Hamilton's 
equations (fTTjl and (fH^ to be 

dU{t)^-pci^{0,t)^l?q{t) , (21) 

The equations of motion (|20|) and H21|l and boundary condition (|19|) are satisfied 
by the bulk solutions, 

/•OO 

(/)(x, t) = {p{p^ - A*^ + ™^) cos{px) - (A(/9^ - + m^) + P'^) s\n{px)) 
Jo 

X {a{p) cos{LUpt) + b{p) siii{LOpt)) , (22) 



qit) = / f3 p {a{p) cos{lj pt) +b{p)sm{iL;pt)) dp , (23) 
Jo 

Where uip ~ \/ m? + p^ and a(/9) and 6(p) are real functions of p. 

These solutions have several interesting features. For large A or /3 the solutions 
for the field become equivalent, after a rescaling of a(p) and 6(/o), to solutions for the 
Dirichlet boundary condition and for small A and /3 the field solutions become the 
same as for a Neumann boundary. These qualitative features of the bulk solutions 
are shared by the solutions for the Robin boundary condition discussed section 1, 
this behaviour is said to 'extrapolate between the Dirichlet and Neumann boundary 
conditions'. The solutions for the dynamic boundary (|22|l and (|23l) contain addi- 
tional interesting features. In the region of the parameter space where p > m there 
exists a value oi p, p = \/ p? ~ vn? at which for any value of /3 and A the solution, 
for this one mode, becomes proportional to sin(pa;), this can be interpreted as a 
resonance effect, at this value of p the field is oscillating with frequency Wp = /i, 
the natural frequency of the boundary oscillator. A second feature in the solu- 
tions not previously observed for just the Robin boundary potential is that when 
\p^ = X{p? — m?) — 0^ the single mode solution becomes proportional to cos(pa;). 
In the next section we will return to these features when we consider the reflection 
cross section of this boundary. 

In addition to the spatially oscillating 'bulk' solutions there exist square inte- 
grable 'boundary bound state' solutions found from the bulk solutions by allowing 
the parameter p to take imaginary values p — ig and requiring that the divergence 
of the field on the left half line arising from the terms proportional to cos]i{gx) and 
smh.(gx) cancel. These solutions have the form 

(j){x, t) = e^'-^ {ug cos{ujgt) + bg sm{ujgt)) , (24) 
q{t) = - ^-^^ {agCos{u)gt) + bgsm{ujgt)) , (25) 

where ujg = ^ w? — g'^, ag and bg are real amplitudes and g \s a positive real 
solution of 

g^ + Xg^ + (p^ - m^)g + X{p^ - m^) -13^=0. (26) 
We can split the parameter space of the theory into regions, using the planes 

• p? ~ m^=0 

• A = 

• / = A(^2 _ ^2) _ ^2 ^ Q 

• J± = |A(//2 _ m^) + ± ^(A2 - 3(Ai2 - m"))^/^ - /J^ = 
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The value of / is the value of the cubic + Xg^ + {p? — m?)g + A(^^ — m?) — f}^ 
at £< = 0, J-i- is the value of the cubic at the local maximum of the function and J_ 
is the value of the cubic at the local minimum. If < 3(/i^ — m^) there are no 
stationary points and J± take complex values. When this occurs we use the real 
part of J± = as the boundary between regions. Figure 2 shows how the parameter 
space of the theory is divided up by these planes. 

It is simple to find how many real solutions of l|26|l exist in each of these regions. 
For fi > m there is a single positive real solution if / < which corresponds to 
regions (a), (b) and (c) in figure 2A, if J+ < 0, i.e. region (c), there are also two 
negative real solutions to H26I) . U p > m and / > then there are no positive real 
solutions, but there is a single real negative solution in regions (d) and (f) and three 
negative solutions in region (e). When fi < m there is a single positive solution 
to if / < 0, i.e. in regions (g), (h), (j) and (k) of figure 2B. In regions (h) 
and (k) there are also two negative real solutions. When /i < m and / > 0, i.e. in 
region (i) there are two positive real solutions and one negative real solution. As 
only positive real solutions of H26|) correspond to classical boundary bound state 
solutions this means there exist one such solution in regions (a), (b), (c), (g), (h), 
(j) and (k). There are no boundary bound states in regions (d), (e) and (f) and 
there are two classical boundary bound state solutions in region (i). 

For H24|l and 1)25(1 to describe time oscillatory states g must be less than m. 
Using (|26|l it is possible to express this requirement as an inequality between the 
parameters of the model, 



If this condition is satisfied the boundary bound states of the system are time 
oscillatory, otherwise they are solutions in which the field diverges as t ^ ±00. 

Classical boundary bound state solutions for the Klein-Gordon field have previ- 
ously been observed for the Robin boundary, as discussed in section 1, where the 
equivalent condition to ((47|l for the bound state to be oscillatory is A < m. It is 
worth noting that the bound state only exists for the attractive Robin boundary, 
i.e. A < 0. If A > then the Robin boundary is repulsive and no bound state 
exists. 

We will now show that the solutions (|22|I - H25() can be written as a superposition 
of orthogonal independent modes of oscillation. Let us define 




(27) 




(28) 



(29) 




(30) 



(31) 



■0(p, X) 



p{p'^ - p^ + m^) cos{px) — {X{p^ — p'^ +m'^ + /3^) sm{px) 



(32) 



xip) 




(33) 
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Figure 2. The various regions of the parameter space of the dy- 
namic boundary with Robin potential theory for (A) fi > m and 
(B) n<m. 
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We can now write the general solutions in terms of the definitions as 

(/)(a;, t) = cl{p)^{p^ ^) cos{ujpt) H x) sm{Ljpt) dp 

Jo \ J 

+ i dgtPg{x)COS{uJgt) + ^'lljg{x)Sm{ujgt)] , (36) 

poo I h( n\ \ 

q{t) ^ J f a{p)x{p) cos(wpt) + —X{p) sm{ujpt) j dp 

+ y2\"'eXQC0s{Ljgt) + —XgSm{ujgt)\ . (37) 
. V '^^ / 

The sum in H36(l and H37|l runs over all positive real roots of H26(l . The functions 
()32|l - ()35|l form an orthonormal set over the half line plus the oscillator, 

■4^{pi,x)i(;{p2,x) dx + x{pi)x{P2) ^S{pi - P2) , (38) 

il;g{x)ip{p,x) dx + XeXip) =^0 , (39) 

g2 (a;) dx -\r XgiXQ2 ~ ^qiQ2 ■ 

(40) 

Relation H38|) holds as p is always positive, H39|l requires that is a root of H26|l 
and uses the relations between roots of a cubic equation^. By substituting 
the general solutions (|36|l and H37|l into the Hamiltonian H14|l and applying the 
othoganality relations l|38|l - (|40|l we can rewrite the Hamiltonian as an infinite sum 
of independent harmonic oscillators each corresponding to one mode of oscillation 
of the field^ 

H^\j^ h{pf + ul~a{pf dp+\Y^bl + ujldl . (41) 

Note that the Hamiltonian l|41|l is bounded below only if uig is real, i.e. if all 
the boundary bound state modes are oscillatory, which occurs only when H27|l is 
satisfied. Thus we must require that the condition H27(l be satisfied to have a 
physical system. 

2.2. Quantum System. By using the definitions (|28|I - H35(I we have described the 
general classical solutions for the field as the linear superposition of independent 
modes of oscillation. The Hamiltonian in this basis can be written as a sum over 
an infinite number of harmonic oscillators each corresponding to one of the modes 



^Let the roots of the cubic + C2 2'^ + ciz + cq be z\, 22 and 23. These roots satisfy the 
relations 212223 = — cq, 2122 + 2123 + 2223 = ci and 21 + 22 + 23 = -C2. 

'^I.e. the integral over all bulk modes plus the sum over the fininte number of boundstate 
modes. 
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of oscillation 141(1 . It is now simple to quantize the system, by analogy to the quan- 
tization of the unbounded field we can define creation and annihilation operators, 



Where there is one pair of creation annihilation operators for each allowed boundary 
bound state. The operators have the usual commutation relations for creation and 
annihilation operators, [a{pi) , a"^ {p2)] = S{pi — P2), [ag,Q;^] = 1, and [a(p),Q;^] = 
[a„at(p)] =0. 

Using these operators we can construct a Fock space, let the vacuum, |0), be the 
state which vanishes when acted on by any of the annihilation operators, Q;(p)|0) = 
0, ag\0) = 0. The one-particle 'bulk' states are created by the action of a^p) on the 
vacuum, Q;^(p)|0} = \p), and one particle boundary states are created by the action 
of a], on the vacuum, q;|,|0) = \g). The one-particle bulk states, |p), represent 
a particle of mass m that at t = —00 is located far to the left of the boundary 
and is travelling with momentum p = p towards the boundary. At time t ~ 00 this 
particle is again located far to the left but is travelling away from the boundary with 
momentum p ~ —p. At some finite time the particle has approached and reflected 
from the boundary. Boundary states, \q), represent particles which are located close 
to the boundary at all times. Annihilation operators act on the one-particle states 
to return the vacuum or zero, a{pi)\p2) — S{pi — P2)|0), agj£'2) = <^eie2|0)- Multi- 
particle states can be constructed from the action of several creation operators on 
the vacuum. We can define the number operator for both the bulk particle states, 
N{p) := a^{p)a{p), and for the boundary states, Ng := a^gUg, from which we can 
construct the renormalized Hamiltonian, 



This is equivalent to the classical Hamiltonian 1(41(1 up to a constant infinite term. 
From equation ((42(1 and the definition of cOp we see that the renormalized Hamilton- 
ian is hermitian if and only if ujg is real. This is the same requirement as that for 
the classical boundary bound states to be oscillatory in time. It has already been 
shown that this condition is equivalent to equation 1(27(1 . Thus provided that the 
condition 1(27(1 between the parameters of the model hold the renormalized Hamil- 
tonian 1(42(1 will be a hermitian operator. If 1(47(1 does not hold the theory has an 
unstable vacuum and cannot be quantized. 

From the general solution for the field ((36(1 and the definition l(12() we find the 
reflection matrix to be 



Note that R{p) is unitary R{p)*R{p) — 1 and that as the parameter p is the 
momentum of incoming particles and as the boundary defines the right most limit 
of the field all incoming particles will have positive momentum, thus it is sufhcient 
for us to consider just these values of p. 

The reflection matrix 1(43(1 has three poles corresponding to the three complex 
roots of 1(26(1 . These poles either occure at positive, purely imaginary, values of 
momentum (p — ig) in which case they correspond to the bound state of the 
system, or they occure at values of momentum in the lower half complex plane. 





(42) 



R{p) 



p (p2 - + m^) - i (A(p2 - + m'^) + P'^) 
p ijP — p? + m?) + i (A(p2 — p? + m?) -t- /S^) 



(43) 
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Figure 3. The total reflection cross sections of (A) the dynamic 
boundary with m = 4, /^ = 5, /3 = 2, A = and (B) the repulsive 
Robin boundary with \ — 



To help interpret these poles we define the total reflection cross section by analogy 
with the optical theorem for bulk scattering, 



where Q{p) = {—i){R{p) — 1) and indicates that we are taking the imagi- 

nary part of <5(p). The total reflection cross section crtot(p), measures the probabil- 
ity that, during the reflection process, the reflecting particle in state \p) is scattered 
into some intermediate state, |^^), which then decays back into the original state*. 
Figure 3A shows the shape of titot for the dynamic boundary in the region of the 
parameter space where p, > m, for the special case when A = 0. For comparison fig- 
ure SB shows the shape of the total reflection cross section of the Robin boundary. 
Both reflection cross sections exhibit a background effect where, for low momentum 
particles, the cross section increases arising from the attraction or repulsion of the 
boundary to the scattering bulk particles. The reflection cross section (|44ll has two 
other important features, the first is a resonance peak around p ~ \/ P^ — fn^ where 
the cross section is large the second is a minimum where the cross section falls to 
zero at Xp^ = \{p? — w?) — (P . These two features correspond to when the classical 
solutions are proportional to sin(jox) and cos(joa;) respectively. Note that Figure 3A 
also exhibits a resonance peak but no minimum, we will discuss this special case in 
a later section. The resonance is only observable in the reflection cross section when 
/i > m as p is positive and real for the scattering particles. Similarly the minimum 
will only appear in the regions of the parameter space where p > m and A < 0, 
i.e. region (a) of figure 2A, or A > and (3^ < A(^^ — m^), regions (d), (e) and (f) 
of figure 2A. For p < m the minimum will appear if A < and > A(/i^ — m^), 
regions (g) and (h) of figure 2B. Figure 4 shows the shape of the reflection cross 
section in the different regions of the parameter space. Note the existence of the 
resonance peak in figures 4A, 4B and 4C and the zero cross section in 4A, after the 
resonance, 4B, before the resonance, and in 4E. 



The total cross section ctot here does not have a direct clear interpretation due to quantum 
mechanical interference between scattering and non-scattering terms, we pursue this interpretation 
through analogy with the bulk case. 




(44) 



p2 (p2 - ^2 _^ ^2)2 _^ f^^p2 _ ^2 _^ ^2) ^ ^2) 
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Figure 4. The total reflection cross section for the boundary os- 
cillator with Robin potential for (A) A = —2, /x = 13, m = 12, 
/3 = 4, (B) A = 2, /X = 13, m = 12, /J = 4, (C) A = 2, = 13, 
m = 12, /? = 8, (D) A = -2, /X = 12, m = 13, /? = 4, (E) A = -2, 
^ = 12, m = 13, /3 = 8, (F) A = 2, /z = 12, m = 13, /3 = 4. 
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The shape of the cross section close to a narrow resonance peak is well approxi- 
mated by the Breit-Wigner formula |^ ^] 

atotip)^(^itOp-toof + \T^^ , (45) 

where uip is the energy of the reflecting particle, is the mass of the resonance 
particle and F is the decay rate of the resonance. By comparing (I45|l and (|44|l close 
to the resonance peak we find that the energy of the resonance particle is equal to fj,, 
the natural frequency of the boundary oscillator. This suggests a new interpretation 
for the resonance peak, if we view the boundary terms in the Hamiltonian, (|14|l . as 
describing a + 1 dimensional field rather than a mechanical system then particles 
of this field will have mass /i. Such particles would be strongly confined to the 
boundary as they are only able to propagate in the domain of the boundary field. 
The resonance peak corresponds to a bulk state, \p) with energy close to the mass, 
/Lt, of the resonance state, \<,), approaching the boundary and, through the linear 
coupling term, P(j)(0,t)q{t), creates a boundary particle in state |<;) which remains 
on the boundary until it decays back into the bulk state \p). 

From 145|l and 144|l it is also possible to calculate an approximate value for F 
which will be valid when the resonance is narrow, 

A 2 • (46) 

As fi approaches m the shape of the resonance peak widens corresponding to an 
increasing decay rate for the boundary particle. The Breit-Wigner formula H45|l is 
only valid for narrow peaks so (|46|) ceases to be a good approximation. 

Resonances associated with the boundary have been previously observed in 

cniEniiTii. 

The zero reflection cross section arises from the combined effects of the boundary 
oscillator and the Robin boundary potential. A similar effect is observed in quantum 
mechanical scattering from a square well potential, see, for example, |521I2S1- Such 
minima in cross sections are often associated with the Ramsauer-Townsend effect 
where scattering of electrons from noble gas atoms show a minimum for certain 
energies of scattering electrons |24| . 

We will now turn our attention to the pole structure of the quantum reflection 
matrix l|43|l . taking each of the regions of the parameter space in turn. In each case 
the reflection matrix has three poles corresponding to the three complex solutions 
of the cubic H26(l . First let us consider the region where ^ > m and A < 0, region 
(a) of figure 2A. Figure 5 shows the analytic structure of the reflection matrix in 
this region, black circles denote poles and white circles denote nodes. In this case 
there in one pole on the positive imaginary axis and two in the lower complex half 
plane. 

Earlier we observed a single classical boundary bound state solution when /i > m 
provided > A(/i^ — m?). As (3"^ and (^^ — m^) are positive and A is negative 
there exists a boundary bound state in this region which corresponds to the pole on 
the positive imaginary axis in figure 5. Earlier in this section we observed that the 
reflection cross section in this region of the parameter space exhibits a resonance 
peak, as illustrated in figure 4A. The unstable resonance state associated with this 
feature can be identified with the pole in the lower right quadrant of figure 5. The 
third pole in the lower left quadrant produces a resonance feature in the cross 
section at p = — y/i^ — however, as mentioned before, p is a positive parameter 
so this resonance is never observed. 

It is illuminating to consider how the poles behave in the weak coupling limit, 
i.e. as /3 — > and A — * simultaneously. In this case the pole corresponding to the 
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Figure 5. The analytic structure of the reflection matrix of the 
dynamic boundary with Robin potential for > m and A < 
(region (a)). 
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Figure 6. The analytic structure of the reflection matrix of the 
dynamic boundary with Robin potential for > m, A > and (A) 
/32 < J_ (region (f)) (B) J_ < /J^ < J+ < A(/x2 - m^), or J_ < 
/32 < X{^^ - to2) < J+ (region (e)), (C) J+ < < X{fj,'' - m^) 
(region (d)). 



bound state approaches p = along the positive imaginary axis. This indicates that 
the state arises from the stationary bulk particle state, the same behaviour as that 
of the bound state of the attractive Robin boundary. Such states can be interpreted 
as particles of the field with energy less than their mass, m. These particles can 
only exist close to an attractive boundary. Thus the boundary in this region of the 
parameter space is attractive. The poles corresponding to the resonance states of 
the system in the weak coupling limit move onto the real p axis at p = p? — w?. 
So when the coupling parameters are reduced the resonance state stabilizes with 
energy equal to the mass of the boundary field, p, as would be expected for a state 
which arises from the particle state of the boundary field. 

We will now consider the region of the parameter space where /i > to, A > 
and 0^ < A(/i^ — m^), regions (d), (e) and (f). Figure 6 illustrates the possible 
configurations of the poles of the reflection matrix in these regions. In each case 
all three poles are located in the lower complex half plane, and at least one pole is 
located on the imaginary axis. 

We know that there are no classical boundary bound state solutions in this 
region of the parameter space, poles located on the negative imaginary axis do not 
represent boundary bound states like those located on the positive imaginary axis. 
We also know that there is a resonance peak in this region. Let us consider the 
behaviour of the poles in the weak coupling limit^, one pole approaches p = along 
the negative imaginary axis, this is the same weak coupling behaviour displayed 
by the pole for the repulsive Robin boundary. Such poles do not carry a physical 



'where A and P are reduced such that < A(/i^ — m^) is always true. 
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Figure 7. The analytic structure of the reflection matrix of the 
dynamic boundary with Robin potential for > m, A > and 
(A)/32 > A(At2 - to2) > J+, or /3'^ > J+ > X{ij,'^ - rn?) (region (b)), 



(B)A(/i^ - m^) <J+ (region (c)). 



interpretation as a valid state of the system. Boundaries whose reflection matrix 
contains poles which have this weak coupling behaviour are repulsive to the particles 
of the field. The remaining poles of the reflection matrix approach p = ±^/Jjfl^^rn? 
in the same way as described in the previous case. As before we can interpret one 
of these poles as being the resonance observed in figure 4B, even though it may be 
located at purely imaginary values of p. Neither of the other poles of the reflection 
matrix have physical interpretations. 

The third case we will consider is for > to, A > and > A(/z^ — m?), regions 
(b) and (c) of figure 2A. Figure 7 shows the possible configurations of the poles of 
R{p) in this region. Again we have a resonance as seen in figure 4C and we also 
have a classical bound state solution, which corresponds to the pole on the positive 
imaginary axis as usual. If we take the low coupling limit, ensuring 0^ > A(/i^ — m^) 
remains valid, we find that this pole approaches p = along the positive imaginary 
axis indicating that the boundary in this region is attractive. We can consider 
how the state behaves if we were to reduce (3^ and hold A constant (or increase 
A and hold f3 constant). In this case the bound state pole passes through p = 
when = A(/x^ — to^). If is decreased further we will pass into a region of the 
parameter space where the boundary is repulsive. We know that this region does 
not have a bound state, thus we conclude that the state corresponding to the bound 
state pole evaporates into a stationary bulk particle state as we make the transition 
between the two regions. The remaining poles behave in the low coupling limit in 
the same manner as described in the previous two cases and one is associated with 
the resonance state. 

The fourth region we will consider is p < m, X < and < A(/u^ — m^), region 
(i) of figure 2B, we know that in this region there arc two classical boundary bound 
state solutions, but there is no resonance state. Figure 8 shows the arrangement 
of poles of R{p) for this region, as expected there are two poles on the positive 
imaginary axis corresponding to the two bound states. 

In the weak coupling limit one of the bound state poles approaches p — along 
the positive imaginary axis, indicating that this state is a boundary bound state in 
the sense of an attractive Robin boundary bound state, it arises from the stationary 
bulk particle state as the coupling is turned on. The presence of this state also 
indicates that the boundary is attractive in this region. The other poles behave 
differently in the weak coupling limit, they move towards the purely imaginary 
values, p = ±y/ p^ — ni^. The energy of these states becomes the same as the mass 
of the boundary field, p. Clearly these states arise from the particle state of the 
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Figure 8. The analytic structure of the reflection matrix of the 
dynamic boundary with Robin potential for /.i < to, A < and 
< A(^^ - TO^) (region (i)). 





o o 



Figure 9. The analytic structure of the reflection matrix of the 
dynamic boundary with Robin potential for /i < m and (A) 0^ > 
J-l_, (regions (g) and (j)), and (B) A(/x^ — m?) < < J+ (regions 
(h) and (k)). 

boundary field in the same manner as the resonances do when /i > m. The pole on 
the positive imaginary axis is interpreted as the stable boundary state arising from 
the boundary field, the other pole lacks a physical interpretation. 

The final parameter space regions we need to consider is when /i < m, and 
0^ > A(/i^ — TO^), which includes everywhere that A > 0, the regions (g), (h), (j) and 
(k) of figure 2B. We know that in these regions there is a single classical boundary 
bound state and no resonance state. Figure 9 shows the possible analytic structure 
of the refiection matrix in these regions. In the weak coupling limit the pole on the 
positive imaginary axis moves to p = \/ y?' — ni?, and carries the interpretation of 
being the state arising from the boundary particle state. A second pole approaches 
p = along the negative imaginary axis, indicating that the boundary is repulsive. 

2.3. The A = special case. In this section we will discuss a special case of the 
system considered previously where A = i.e. where the Robin boundary potential 
is removed leaving the field just coupled to the oscillator at the boundary. 

Figure 10 shows the possible real solutions of (|26|l for A = and for different 
values of m, /i and /?, clearly there is one and only one real positive solution for 
any allowed choice of the parameters. 

The condition on the parameters of the theory required for the bound state 
corresponding to this solution to be time oscillatory (and for the Hamiltonian to 
be bounded below) l|77|l becomes 




(47) 
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Figure 10. Graphs of the cubic + — m^) — when (A) 
/X > m, (B) ^ < m and > ^(3(to^ - A^^))^^^ and (C) fi < m 
and /32 < ^(3(m2-^2))3/2_ 
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Figure 1 1 . The analytic structure of the reflection matrix of the 
dynamic boundary for /i > m, X — 0. 



We note that when m = equation (|47|l can only be satisfied by /3 = 0, i.e. 
a massless field in the bulk can only be coupled to a boundary oscillator in the 
presence of a (repulsive) Robin boundary potential. 

In the quantized system the reflection matrix 1431) simplifles to 

RiP) = f.'^l^"'2'"C (48) 

The reflection matrix (|48|) has three poles, corresponding to the three complex roots 
of H2t)|l . one of these poles always occurs at positive, purely imaginary, values of 
momentum (p = ig) and corresponds to the classical bound state of the system, the 
features of which have been discussed in previously. The remaining two poles occur 
at values of momentum in the lower half complex plane. To help interpret these 
poles we will initially assume that /i > to, in which case the poles are arranged as 
shown in Figure 11. The dotted circle describes the locus of points with the same 
modulus as the location of the bound state pole. 

Figure 3A shows the shape of CTtot for the dynamic boundary in the region of the 
parameter space where /i > to and A = 0. Figure 3A exhibits a resonance peak, the 
pole in the lower right quadrant of figure 11 corresponds to the resonance state. 

We now have a full description of two of the poles of R{p) in the region where 
H > m. The third pole of R{p) carries a similar interpretation to that just presented 
for the resonance. If we were to continue figure 3A into negative values for p 
we would observe a second resonance peak corresponding to the third pole of the 
reflection matrix. However, as mentioned before this region of the reflection matrix, 
and thus the cross section is never probed. 

Before we conclude our discussion of this region of the parameter space let us 
consider the behaviour of the poles of the reflection matrix in the limit /3 ^ 0, 
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Figure 12. The analytic structure of the reflection matrix for the 
dynamic boundary for (A) fi < m and > ^(3(m^ — fP))^/^ and 



(B) fi < m and < 



!))3/2 



where the osciflator decouples from the field. From H48|l we see that in this limit 
the pole corresponding to the bound state approaches g ^ from above. The 
energy of this state in the /9 — > limit is equal to the mass of the bulk particle, 
TO. This is the same behaviour as for the boundary bound state of the attractive 
Robin boundary potential in the limit A — > and indicates the bound state arises 
from the stationary bulk particle state. The ability for a boundary to support such 
a state indicates it is attractive to the particles of the field. In the limit /3 — > the 
poles associated with the resonances stabilise and have energy equal to the mass of 
the particles of the boundary oscillator, /it. This further supports the interpretation 
of the resonance particles as particles of the boundary field. 

Let us now consider the region of the parameter space where fj, < m. Figure 12 
show the arrangement of poles and nodes of the reflection matrix in this region. 

The reflection cross section can be calculated and plotted as previously, however 
no additional features beyond background effects at low momentum are observed. 
To interpret these poles let us observe their behaviour as /9 — > 0. When < to 
we find that the bound state pole, being the one located at p — ig on the positive 
imaginary axis, no longer goes to g = as the boundary coupling parameter is 
reduced. Instead its limiting energy is /i, the mass of the boundary particles. As 
in the previous case there is also a second pole which also has an energy /j, in the 
weak coupling limit, but as it is located on the negative imaginary axis it carries no 
physical meaning. The final pole approaches g = from below as P is reduced, this 
is the same behaviour as the pole in the reflection matrix of the repulsive Robin 
boundary. In both cases the pole does not correspond to any kind of bound state. 

Clearly the spectrum of boundary states of the dynamic boundary is very dif- 
ferent in these two regions of the parameter space. For /i > m there is a resonance 
state corresponding to the particles of the boundary field and a Robin type bound- 
ary bound state arising from the attraction of the bulk particles to the boundary. 
However, when /i is reduced past m the Robin type bound state vanishes, just as for 
the Robin boundary when A is increased past zero. The resonance state becomes a 
new, stable, boundary bound state, classically this state can be thought of as the 
oscillator driving the field at the boundary but at a frequency too low to allow the 
energy to be dissipated away. The absence of a Robin type boundary bound state 
suggests that the boundary oscillator in this region is repulsive to the particles of 
the field. 
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3. Conclusion 

In this paper we considered the massive Klein-Gordon field coupled to a boundary 
oscillator, with and without an additional Robin boundary potential. We showed 
that these systems can be solved classically and that the classical solutions can 
be written as a superposition of independent modes of oscillation. We observed 
classical boundary bound state solutions in some regions of the parameter space, 
the requirement that these solutions be time oscillatory is seen to be the same as 
the requirement that the Hamiltonian is bounded below. We develop a condition 
relating the parameters of the theories which, when satisfied, ensures that bound 
state is oscillatory and thus that the system is physical. 

Writing the classical solutions as superposition of independent modes allows the 
system to be quantized using canonical methods. The quantum reflection matri- 
ces for the boundaries are found from the two point function of the field. We 
observed several interesting features in the total cross section of the refiection pro- 
cess, including resonances and Ramsauer-Townsend effects. Although resonances 
associated with boundaries have previously been observed P El HUJ they have 
only previously been associated with boundaries containing additional degrees of 
freedom in jl4| . Ramsauer-Townsend effects have not previously been observed in 
boundary field theories. 

Quantum boundary bound states were observed in several regions of the pa- 
rameter space. Some of these bound states were seen to arise from the states of 
the boundary oscillator, this behaviour has not previously been observed. Other 
boundary bound states arise in a similar manner to those of the attractive Robin 
boundary. 

Similar systems to the ones described in this paper have previously been investi- 
gated. The coupling of one or more harmonic oscillators at fixed points in the bulk 
of a 1+1 dimensional, massless, classical Klein-Gordon field is discussed in |10| . 
and the coupling of an harmonic oscillator to a massless scalar field in a spherical 
reflecting cavity is discussed in [111 IT^ where it is used to model the behaviour of 
an excited atom in a cavity. In the massless scalar fleld coupled to a boundary 
oscillator is studied as a model of brane-bulk interaction in a braneworld universe. 
There also exist several papers dealing with quantum mechanical systems linked 
the some dissipative medium, based on the appoach of ,25^ where the dissapative 
medium is modeled by coupling the system to an inflnite number of harmonic os- 
cillators. This modelling of the dispersive medium should be equivalent to the 
coupling of the system to the Klein-Gordon field, which contains an infinite num- 
ber of harmonic oscillators, however, a correspondence between the results of this 
paper and those from dissapative quantum mechanics has yet to be found. For a 
recent treatment of the harmonic oscillator using the methods of see |26| . 

It seems that similar results to those presented in this paper should be obtainable 
for fields bounded by arbitrary linear mechanical systems. Of particular interest 
would be the existence and energy of boundary bound states and resonances for 
such systems. Other possibly interesting extensions include coupling two or more 
non-interacting bulk fields through a boundary or reformulating the present system 
to describe a point interaction on the whole line, or in a volume. In the latter case 
we would expect to make contact with the results of the earlier papers mentioned 
above. 



20 



KLEIN-GORDON FIELD WITH BOUNDARY OSCILLATOR 



Acknowledgments 

The author would hke to thank Gustav Dehus, Zoltan Bajnok, Laszlo PaUa, 
Gabor Takacs, Evgeni Sklyanin and Bernard Kay for useful and interesting discus- 
sions. The author was supported by a PPARC studentship whilst undertaking this 
research. 



References 

[1] S. Ghoshal and A. B. Zamolodchikov, Boundary S Matrix and Boundary State in Two- 
Dimensional Integrable Field Theory, Int. J. Mod. Phys. A 9 (1994) 3841, Erratum-ibid 9 
(1994) 4353. preprint hep-th/9306002 

[2] E. K. Sklyanin, Boundary Conditions for Integrable Quantum Systems, J. Phys. A 21 (1988) 
2375. 

[3] Z. Bajnok, G. Bolim and G. Takacs, On Perturbative Quantum Field Theory with Boundary, 
Nucl. Phys. B 682 (2004) 585, preprint hep-th/0309119 

[4] V. Bazhanov, S. Lukyanov and A. Zamolodchikov, On Nonequalibrium States in QFT Model 
with Boundary Interaction, Nucl. Phys. B 546 (1999) 529, preprint hep-th/9812091 

[5] P. Baseilhac and G. W. Delius, Coupling Integrable Field Theories to Mechanical Systems at 
the Boundary, J. Phys. A 34 (2001) 8259, preprint hep-th/0106275 

[6] P. Baseilhac, G. W. Delius, A. George, Coupling the Sine-Gordon Field Theory to a Mechan- 
ical System at the Boundary, in Statistical Field Theories, ed A. Cappelli and G. Mussardo 
(Klewer 2002) preprint nlin.SI/0201007 

[7] P. Baseilhac and K. Kuizumi, Sine-Gordon Field Theory on the Half-Line with Quantum 
Boundary Degrees of Freedom, Nucl. Phys. B 649 (2003) 491, preprint hep-th/0208005 

[8] R. Napomechie, The Boundary Supersymmetric Sine-Gordon Model Revisited, Phys. Lett. B 
509 (2001) 183, preprint hep-th/0103029 

[9] M. Fuentes, A. Lopez, E. Fradkin and E. Moreno, Bozonization Rules in 1/2 -|- 1 Dimensions, 
Nucl. Phys. B 450 (1995) 603, preprint cond-mat/9502076 
[10] S. A. Chorosavin, Simple Exactly Solvable Models Based on Finite Rank Perturbation Meth- 
ods I: D'Alembert-Kirchhoff-like Formulae, preprint maths. DS/0301137. 

S. A. Chorosavin, Simple Exactly Solvable Models Based on Finite Rank Perturbation Meth- 
ods II: Resolvents Formulae, preprint maths. DS/0301137. 

S. A. Chorosavin, Simple Exactly Solvable Models Based on Finite Rank Perturbation Meth- 
ods III: Linear Friction as Radiation Reaction, preprint math-ph/0307029 
S. A. Chorosavin, Simple Exactly Solvable Models of One-Dimensional Scalar Fields with 
Concentrate Factors, preprint math-ph/0311045 
[11] N. P. Andion, A. P. C. Malbouisson and A. Mattos Neto, An Exact Approach to the Oscil- 
lation Radiation Problem in an Arbitrary Large Cavity, J. Phys. A 34 (2001) 3735, preprint 
physics/0001009 

[12] G. Flores-Hidalgo, A. P. C. Malbouisson and Y. W. Milla, Stability of Excited Atoms in 
Small Cavities, Phys. Rev. A 65 (2002) 063414. 

[13] G. Flores-Hidalgo and A. P. C. Malbouisson, Time Evolution of Quantum Systems in Micro- 
cavities and in Free Space - a Non-Perturbative Approach, preprint hep-th/0308198 

[14] P. Binetruy, M. Bucher and C. Carvalho, Models of the Brane-Bulk Interaction: Toward Un- 
derstanding Braneworld Cosmological Perturbation, Phys. Rev. D 70 (2004) 043509, preprint 
hep-th/0403154 

[15] S. A. Fulling, Aspects of Quantum Field Theory in Curved Spacetimes (Cambridge University 
Press, 1989). 

[16] Z. Bajnok, G. Bohm and G. Takacs, Boundary Reduction Formula, J. Phys. A 35 (2002) 

9333, preprint hep-th/0207079 
[17] M. E. Peskin and D. V. Schroder, An Introduction to Quantum Field Theory (Addison- Wesley 

1995). 

[18] S. Weinberg, The Quantum Theory of Fields (Cambridge University Press, 1995). 

[19] P. Mosconi, G. Mussardo and V. Riva, Boundary Quantum Field Theories with Infinite 

Resonance States, Nucl. Phys. B 621 (2002) 571, preprint hep-th/0107082 
[20] Z. Bajnok, L. Palla and G. Takacs, The Spectrum of Boundary Sine-Gordon Theory, preprint 

hep-th/0211132 

[21] G. Breit and E. P. Wigner, Capture of Slow Neutrons, Phys. Rev. 49 (1936) 519. 

[22] A. I. M. Rae, Quantum Mechanics 3rd Ed. (loP 1992). 

[23] R. C. Greenhow, Introductory Quantum Mechanics (loP 1990). 

[24] J. Holtsmark, Der Ramsauereffekt im Argon, Z. Phys. 55 (1929) 437. 



KLEIN-GORDON FIELD WITH BOUNDARY OSCILLATOR 



21 



[25] A. O. Caldeira and A. J. Loggctt, Influence of Dissipation on Quantum Tunneling in Macro- 
scopic Systems, Phys. Rev. Lett. 46 (1981) 211. 

A. O. Caldeira and A.J. Leggett, Quantum Tunneling in Dissipative Systems, Ann. Phys. 
149 (1983) 374. 

[26] R. Rosenfelder, Structure Function of a Damped Harmonic Oscillator, Phys. Rev. C 68 
(2003) 034602, preprint nucl-th/0303003, 



